Microsatellites used in forensics are in regions enriched for trait-associated variants.

The 20 short tandem repeat (STR) loci of the combined DNA index system (CODIS) are the basis of the vast majority of forensic genetics in the United States. One argument for permissive rules about the collection of CODIS genotypes is that the CODIS loci are thought to contain little information about ancestry or traits. However, in the past 20 years, a growing field has identified hundreds of thousands of genotype-trait associations. Here, we conduct a survey of the landscape of such associations surrounding the CODIS loci as compared with non-CODIS STRs. Although this study cannot establish or quantify associations between CODIS genotypes and phenotypes, we find that the regions around the CODIS loci are enriched for both known pathogenic variants (> 90th percentile) and for trait-associated SNPs identified in genome-wide association studies (GWAS) (≥ 95th percentile in 10kb and 100kb flanking regions), compared with other random sets of autosomal tetranucleotide-repeat STRs.

Improving biology faculty diversity through a co-hiring policy and faculty agents of change.

Persons Excluded due to Ethnicity and Race (PEERs) remain underrepresented in university faculties, particularly in science, technology, engineering, math and medicine (STEMM) fields, despite increasing representation among students, and mounting evidence supporting the importance of PEER faculty in positively impacting both scientific and educational outcomes. In fact, the ratio of PEER faculty to students has been steadily dropping since 2000. In our case study, we examine the factors that explain creation of an unusually diverse faculty within a biology department. We analyzed nearly 40 years of hiring data in the study department and show that this department (the study department), historically and currently, maintains a significantly higher proportion of PEERs on faculty as compared to two national datasets. Additionally, we identify factors that contributed to hiring of PEERs into tenure and tenure-track positions. We observed a significant increase in the hiring of PEERs concurrent with the implementation of a co-hiring policy (p = 0.04) which allowed a single search to make two hires when at least one candidate was a PEER. In contrast, three key informants at sister departments reported that co-hiring policies did not result in PEER hires, but instead different practices were effective. In line with one of these practices, we observe a possible association between search committees with at least one PEER member and PEER hiring (p = 0.055). Further, the presence of particular faculty members (Agents of Change) on search committees is associated with PEER hiring. In this case study the combination of a co-hire policy based on the principle of interest-convergence to redress hiring inequities, along with the presence of agents of change, increased faculty PEER representation in STEMM departments.

Microsatellites used in forensics are located in regions unusually rich in trait-associated variants.

The 20 short tandem repeat (STR) markers of the combined DNA index system (CODIS) are the basis of the vast majority of forensic genetics in the United States. One argument for permissive rules about the collection of CODIS genotypes is that the CODIS markers are thought to contain information relevant to identification only (such as a human fingerprint would), with little information about ancestry or traits. However, in the past 20 years, a quickly growing field has identified hundreds of thousands of genotype-trait associations. Here we conduct a survey of the landscape of such associations surrounding the CODIS loci as compared with non-CODIS STRs. We find that the regions around the CODIS markers are enriched for both known pathogenic variants (>90th percentile) and for SNPs identified as trait-associated in genome-wide association studies (GWAS) (≥95th percentile in 10kb and 100kb flanking regions), compared with other random sets of autosomal tetranucleotide-repeat STRs. Although it is not obvious how much phenotypic information CODIS would need to convey to strain the "DNA fingerprint" analogy, the CODIS markers, considered as a set, are in regions unusually dense with variants with known phenotypic associations.

Analyses of Genome Regulatory Evolution Following Whole-Genome Duplication Using the Phylogenetic EVE Model.

Whole-genome duplications (WGDs) are important in shaping the evolution of complex genomes, including rewiring of genome regulation. To address key questions about how WGDs impact the evolution of genome regulation, we need to understand the relative importance of selection versus drift and temporal evolutionary dynamics. One promising class of statistical models that can help address such questions are phylogenetic Ornstein-Uhlenbeck (OU) models.Here we present a computational pipeline for the comparative phylogenetic analyses of genome regulation using an OU model. We have implemented this model in R and provide a step-by-step protocol for the use of this model, including example scripts and simulated test data. We provide the nonspecialist a brief overview of how this model works and how to perform tests for signatures of selection on genome regulation as well as power simulations to aid in experimental design and interpretation of results. We believe that these resources could help polyploidy research move forward in an era of rapidly increasing functional genomics data across the tree of life.

Associations between forensic loci and expression levels of neighboring genes may compromise medical privacy.

A set of 20 short tandem repeats (STRs) is used by the US criminal justice system to identify suspects and to maintain a database of genetic profiles for individuals who have been previously convicted or arrested. Some of these STRs were identified in the 1990s, with a preference for markers in putative gene deserts to avoid forensic profiles revealing protected medical information. We revisit that assumption, investigating whether forensic genetic profiles reveal information about gene-expression variation or potential medical information. We find six significant correlations (false discovery rate = 0.23) between the forensic STRs and the expression levels of neighboring genes in lymphoblastoid cell lines. We explore possible mechanisms for these associations, showing evidence compatible with forensic STRs causing expression variation or being in linkage disequilibrium with a causal locus in three cases and weaker or potentially spurious associations in the other three cases. Together, these results suggest that forensic genetic loci may reveal expression levels and, perhaps, medical information.

Elucidating gene expression adaptation of phylogenetically divergent coral holobionts under heat stress.

As coral reefs struggle to survive under climate change, it is crucial to know whether they have the capacity to withstand changing conditions, particularly increasing seawater temperatures. Thermal tolerance requires the integrative response of the different components of the coral holobiont (coral host, algal photosymbiont, and associated microbiome). Here, using a controlled thermal stress experiment across three divergent Caribbean coral species, we attempt to dissect holobiont member metatranscriptome responses from coral taxa with different sensitivities to heat stress and use phylogenetic ANOVA to study the evolution of gene expression adaptation. We show that coral response to heat stress is a complex trait derived from multiple interactions among holobiont members. We identify host and photosymbiont genes that exhibit lineage-specific expression level adaptation and uncover potential roles for bacterial associates in supplementing the metabolic needs of the coral-photosymbiont duo during heat stress. Our results stress the importance of integrative and comparative approaches across a wide range of species to better understand coral survival under the predicted rise in sea surface temperatures.

Comparative regulomics supports pervasive selection on gene dosage following whole genome duplication.

BACKGROUND: Whole genome duplication (WGD) events have played a major role in eukaryotic genome evolution, but the consequence of these extreme events in adaptive genome evolution is still not well understood. To address this knowledge gap, we used a comparative phylogenetic model and transcriptomic data from seven species to infer selection on gene expression in duplicated genes (ohnologs) following the salmonid WGD 80-100 million years ago. RESULTS: We find rare cases of tissue-specific expression evolution but pervasive expression evolution affecting many tissues, reflecting strong selection on maintenance of genome stability following genome doubling. Ohnolog expression levels have evolved mostly asymmetrically, by diverting one ohnolog copy down a path towards lower expression and possible pseudogenization. Loss of expression in one ohnolog is significantly associated with transposable element insertions in promoters and likely driven by selection on gene dosage including selection on stoichiometric balance. We also find symmetric expression shifts, and these are associated with genes under strong evolutionary constraints such as ribosome subunit genes. This possibly reflects selection operating to achieve a gene dose reduction while avoiding accumulation of "toxic mutations". Mechanistically, ohnolog regulatory divergence is dictated by the number of bound transcription factors in promoters, with transposable elements being one likely source of novel binding sites driving tissue-specific gains in expression. CONCLUSIONS: Our results imply pervasive adaptive expression evolution following WGD to overcome the immediate challenges posed by genome doubling and to exploit the long-term genetic opportunities for novel phenotype evolution.

Illuminating Women's Hidden Contribution to Historical Theoretical Population Genetics.

While productivity in academia is measured through authorship, not all scientific contributors have been recognized as authors. We consider nonauthor "acknowledged programmers" (APs), who developed, ran, and sometimes analyzed the results of computer programs. We identified APs in Theoretical Population Biology articles published between 1970 and 1990, finding that APs were disproportionately women (P = 4.0 × 10-10). We note recurrent APs who contributed to several highly-cited manuscripts. The occurrence of APs decreased over time, corresponding to the masculinization of computer programming and the shift of programming responsibilities to individuals credited as authors. We conclude that, while previously overlooked, historically, women have made substantial contributions to computational biology. For a video of this abstract, see: https://vimeo.com/313424402.

Fitting the Balding-Nichols model to forensic databases.

Large forensic databases provide an opportunity to compare observed empirical rates of genotype matching with those expected under forensic genetic models. A number of researchers have taken advantage of this opportunity to validate some forensic genetic approaches, particularly to ensure that estimated rates of genotype matching between unrelated individuals are indeed slight overestimates of those observed. However, these studies have also revealed systematic error trends in genotype probability estimates. In this analysis, we investigate these error trends and show how they result from inappropriate implementation of the Balding-Nichols model in the context of database-wide matching. Specifically, we show that in addition to accounting for increased allelic matching between individuals with recent shared ancestry, studies must account for relatively decreased allelic matching between individuals with more ancient shared ancestry.

Phylogenetic ANOVA: The Expression Variance and Evolution Model for Quantitative Trait Evolution.

A number of methods have been developed for modeling the evolution of a quantitative trait on a phylogeny. These methods have received renewed interest in the context of genome-wide studies of gene expression, in which the expression levels of many genes can be modeled as quantitative traits. We here develop a new method for joint analyses of quantitative traits within- and between species, the Expression Variance and Evolution (EVE) model. The model parameterizes the ratio of population to evolutionary expression variance, facilitating a wide variety of analyses, including a test for lineage-specific shifts in expression level, and a phylogenetic ANOVA that can detect genes with increased or decreased ratios of expression divergence to diversity, analogous to the famous Hudson Kreitman Aguadé (HKA) test used to detect selection at the DNA level. We use simulations to explore the properties of these tests under a variety of circumstances and show that the phylogenetic ANOVA is more accurate than the standard ANOVA (no accounting for phylogeny) sometimes used in transcriptomics. We then apply the EVE model to a mammalian phylogeny of 15 species typed for expression levels in liver tissue. We identify genes with high expression divergence between species as candidates for expression level adaptation, and genes with high expression diversity within species as candidates for expression level conservation and/or plasticity. Using the test for lineage-specific expression shifts, we identify several candidate genes for expression level adaptation on the catarrhine and human lineages, including genes putatively related to dietary changes in humans. We compare these results to those reported previously using a model which ignores expression variance within species, uncovering important differences in performance. We demonstrate the necessity for a phylogenetic model in comparative expression studies and show the utility of the EVE model to detect expression divergence, diversity, and branch-specific shifts.

Modeling gene expression evolution with an extended Ornstein-Uhlenbeck process accounting for within-species variation.

Much of the phenotypic variation observed between even closely related species may be driven by differences in gene expression levels. The current availability of reliable techniques like RNA-Seq, which can quantify expression levels across species, has enabled comparative studies. Ornstein-Uhlenbeck (OU) processes have been proposed to model gene expression evolution as they model both random drift and stabilizing selection and can be extended to model changes in selection regimes. The OU models provide a statistical framework that allows comparisons of specific hypotheses of selective regimes, including random drift, constrained drift, and expression level shifts. In this way, inferences may be made about the mode of selection acting on the expression level of a gene. We augment this model to include within-species expression variance, allowing for modeling of nonevolutionary expression variance that could be caused by individual genetic, environmental, or technical variation. Through simulations, we explore the reliability of parameter estimates and the extent to which different selective regimes can be distinguished using phylogenies of varying size using both the typical OU model and our extended model. We find that if individual variation is not accounted for, nonevolutionary expression variation is often mistaken for strong stabilizing selection. The methods presented in this article are increasingly relevant as comparative expression data becomes more available and researchers turn to expression as a primary evolving phenotype.

The influence of relatives on the efficiency and error rate of familial searching.

We investigate the consequences of adopting the criteria used by the state of California, as described by Myers et al. (2011), for conducting familial searches. We carried out a simulation study of randomly generated profiles of related and unrelated individuals with 13-locus CODIS genotypes and YFiler® Y-chromosome haplotypes, on which the Myers protocol for relative identification was carried out. For Y-chromosome sharing first degree relatives, the Myers protocol has a high probability (80~99%) of identifying their relationship. For unrelated individuals, there is a low probability that an unrelated person in the database will be identified as a first-degree relative. For more distant Y-haplotype sharing relatives (half-siblings, first cousins, half-first cousins or second cousins) there is a substantial probability that the more distant relative will be incorrectly identified as a first-degree relative. For example, there is a 3~18% probability that a first cousin will be identified as a full sibling, with the probability depending on the population background. Although the California familial search policy is likely to identify a first degree relative if his profile is in the database, and it poses little risk of falsely identifying an unrelated individual in a database as a first-degree relative, there is a substantial risk of falsely identifying a more distant Y-haplotype sharing relative in the database as a first-degree relative, with the consequence that their immediate family may become the target for further investigation. This risk falls disproportionately on those ethnic groups that are currently overrepresented in state and federal databases.

Familial identification: population structure and relationship distinguishability.

With the expansion of offender/arrestee DNA profile databases, genetic forensic identification has become commonplace in the United States criminal justice system. Implementation of familial searching has been proposed to extend forensic identification to family members of individuals with profiles in offender/arrestee DNA databases. In familial searching, a partial genetic profile match between a database entrant and a crime scene sample is used to implicate genetic relatives of the database entrant as potential sources of the crime scene sample. In addition to concerns regarding civil liberties, familial searching poses unanswered statistical questions. In this study, we define confidence intervals on estimated likelihood ratios for familial identification. Using these confidence intervals, we consider familial searching in a structured population. We show that relatives and unrelated individuals from population samples with lower gene diversity over the loci considered are less distinguishable. We also consider cases where the most appropriate population sample for individuals considered is unknown. We find that as a less appropriate population sample, and thus allele frequency distribution, is assumed, relatives and unrelated individuals become more difficult to distinguish. In addition, we show that relationship distinguishability increases with the number of markers considered, but decreases for more distant genetic familial relationships. All of these results indicate that caution is warranted in the application of familial searching in structured populations, such as in the United States.

Detecting coevolution through allelic association between physically unlinked loci.

Coevolving interacting genes undergo complementary mutations to maintain their interaction. Distinct combinations of alleles in coevolving genes interact differently, conferring varying degrees of fitness. If this fitness differential is adequately large, the resulting selection for allele matching could maintain allelic association, even between physically unlinked loci. Allelic association is often observed in a population with the use of gametic linkage disequilibrium. However, because the coevolving genes are not necessarily in physical linkage, this is not an appropriate measure of coevolution-induced allelic association. Instead, we propose using both composite linkage disequilibrium (CLD) and a measure of association between genotypes, which we call genotype association (GA). Using a simple selective model, we simulated loci and calculated power for tests of CLD and GA, showing that the tests can detect the allelic association expected under realistic selective pressure. We apply CLD and GA tests to the polymorphic, physically unlinked, and putatively coevolving human gamete-recognition genes ZP3 and ZP3R. We observe unusual allelic association, not attributable to population structure, between ZP3 and ZP3R. This study shows that selection for allele matching can drive allelic association between unlinked loci in a contemporary human population, and that selection can be detected with the use of CLD and GA tests. The observation of this selection is surprising, but reasonable in the highly selected system of fertilization. If confirmed, this sort of selection provides an exception to the paradigm of chromosomal independent assortment.

One-stage design is empirically more powerful than two-stage design for family-based genome-wide association studies.

Finding a genetic marker associated with a trait is a classic problem in human genetics. Recently, two-stage approaches have gained popularity in marker-trait association studies, in part because researchers hope to reduce the multiple testing problem by testing fewer markers in the final stage. We compared one two-stage family-based approach to an analogous single-stage method, calculating the empirical type I error rates and power for both methods using fully simulated data sets modeled on nuclear families with rheumatoid arthritis, and data sets of real single-nucleotide polymorphism genotypes from Centre d'Etude du Polymorphisme Humain pedigrees with simulated traits. In these analyses performed in the absence of population stratification, the single-stage method was consistently more powerful than the two-stage method for a given type I error rate. To explore the sources of this difference, we performed a case study comparing the individual steps of two-stage designs, the two-stage design itself, and the analogous one-stage design.